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Abstract 

We investigate the time evolution of a quark antiquark plasma by solving numerically the 
relativistic transport equations derived on the Hartree level from the Nambu-Jona-Lasinio model. 
We find that the phase transition in the expanding quark antiquark plasma is different as compared 
to that in a static plasma. The expansion competes with the transition and finally quark droplets 
will be formed which subsequently hadronizes. These findings raise the question whether static 
thermal models can make at all any prediction about signals of that transition . 
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The search for the Quark Gluon Plasma (QGP) is one of the challenges of present day nuclear 
physics. It is the driving motivation of the lead beam experiments which started this year at CERN. 
There is first of all the question whether this phase transition can indeed be achieved with accelerated 
nuclei. The even bigger challenge is, however, to investigate which signals bear witness to the existence 
of the plasma, i.e., cannot be created by the dense hadronic medium which is present after the QGP 
has ceased to exist and the quarks are again bound in the hadrons which will finally be detected. 

Up to now almost all investigations have tacitly assumed that the plasma is contained in a box 
whose boundaries can be moved that slowly that the systems is thermally equilibrated. Then, knowing 
the equation of state, the thermodynamical variables are sufficient to describe the transition to the 
hadronic world (for a review of the thermal scenario we refer to [[!]]). The reality, however, is most 
probably quite opposite. We have a gas of almost massless quarks at a temperature of more then 100 
MeV and nothing which confines the plasma to a certain space region. Hence the system expands and 
we are confronted with the not yet investigated question how the phase transition takes place in the 
expanding non thermal system. 

This question cannot be addressed in the framework of thermodynamics but requires a transport 
theory. In principle a transport theory may be derived from each Lagrangian and consequently also 
from the QCD Lagrangian. The attempts which have been made, however, do not come even close 
to an equation which can be used for numerical investigations |J. Therefore one has to rely on phe- 
nomenological Lagrangians. One of these is the Nambu-Jona-Lasinio Lagrangian (NJL)||. In this 
Lagrangian all gluonic degrees of freedom are integrated out. The interaction between the quarks, 
however, is quite reliable modeled, giving the masses of the known hadrons with an acceptable pre- 
cision. Its drawback is that it has no confinement, i.e., cannot describe properly the hadronization. 
Despite of this fact, this Lagrangian seems to be appropriate to start with the investigation of the 
question how a phase transition takes place in an expanding plasma. The results can be compared 
with that of the string models which recently also conjectured that mini plasmas have to be formed 
during the expansion of the system |Q] . There is no doubt that this approach can be regarded only as 
a first step towards the understanding of the expansion of a quark gluon plasma. 

By using the loop expansion approach of the two-point source connected generating functional p] of 
the NJL Lagrangian, Zhang and Wilets || derived a microscopic transport theory to describe the chiral 
dynamics in high energy heavy-ion collisions. This transport theory which describes the time evolution 
of the Green's function is based on the the closed time-path Green's functions formalism At this 
level it seems, however, not yet possible to perform calculations. If one assumes that the dissipation 
parts of the Green's functions are small compared to the dispersive parts (the so-called quasiparticle 
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limit), the quasiparticle energy spectra of quarks and mesons are obtained from the equations giving 
the poles of the Green's functions, i.e., the zeros of the dispersive parts. Consequently, the Green's 
functions representation can be converted into phase space densities 

In this work we focus on the transport equation without collision terms. For the NJL model, the 
Vlasov equation for the quark density distribution /(r, p,t) in the Hartree level reads as follow (6|, |j 
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where E is the relativistic quark energy y/p 2 + M| and M c is the constituent quark mass which is 
determined by the gap equation ||] : 

M c (r, t)=m + g 2 M c (r, t) T - A ^4f ~ /(r, P, t) - f(r, p, t) ) . (2) 

J o VP + M*{r,t) [2,Tiy 

In Eq.(^) to is the quark mass (to = m u = rrid = 4 MeV), N c = 3 is the color number, Nf = 2 is 
the flavor number, g is the NJL four-point interaction coupling constant |J and A is the cutoff in 
momenta. In order to reproduce the pion mass and the pion decay constant the parameters g and A 
are fixed to 0.48 fm and 820 MeV respectively H . 

The solution of Eq.(|l]) is not trivial due to the dependence of the energy E on the distribution / via 
the constituent quark mass (c.f., Eq.(|2|)). For this reason we have to use some approximation scheme. 
In this work we use the test particles method which has been shown to be powerful in intermediate 
energy heavy ion collisions p. This approach consists to approximate the distribution function / (/) 
by a crowd of "numerical" quarks (anti-quarks). Namely, 

A 

/(r, p, t) = w ]T 5 s (r - T t {t)) 5 3 (p - Pi (t)) 

l A (3) 
/>, p, t) = w ^ 3 (r - ii{t)) 5 3 (P ~ Pi(t)) 

where A and A are the number of test-particles and anti-particles respectively, while w is a normal- 
ization factor which is related to the number B of the real quarks of system by 

J d 3 rd 3 p (/(r, p, t) - f(r, p, t)) = w(A — A) = B . (4) 

To satisfy Eq.(|]), the position and momentum of the i-th. numerical quark (anti-quark) should obey 
the relativistic Hamilton's equations. Namely, 



i l = p i /E i , pi = -V r Ei , Ei = <JpKt) + M*(Ti,t) 

(5) 



r; = pi/Ei , ^ = -V r Ei , Ei = Jp 2 (t) + M 2 (t t ,t) 
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where dots indicate time differentiations. We have investigated numerically the equations of motion 
(||) together with the gap equation (|2|) in a self-consistent manner. We have performed our calculations 
on a 3 dimensions 41 point lattice with a a = 0.25 fm mesh size. The initial conditions are taken as 
follows : we have considered A test-quarks and A test anti-quarks (see below for the numerical values) 
distributed randomly within a sphere of radius tq = 1.05 fm. The position of each test particle is 
assigned by choosing: 

r = r (xi) 1/3 ; cos(0) = 1 - 2x 2 ; 4> = 2vnr 3 (6) 

where x%,X2 and X3 are random numbers between and 1. The Cartesian coordinates of the test 
particle are rcos^sin#, r sin (j) sin 6 and rcos9. The momentum of the test particles are choosen 
randomly according to a Fermi distribution. It means that initially the system is in equilibrium. The 
coordinates of each test-particle momentum are \p x = p cos <j) p sin 6 p , p y = p sin cf> p sin 9 P , p z = p cos 9 P ] 
where the angles are choosen randomly as above while p is determined from 

P E /— fj, /"A E / — [± 

p 2 dp' (1 + e^T^y 1 I / p' 2 dp' (1 + e^)" 1 = x x (7) 
J 



In the last equation x\ is a random number, E p > = yp" 2 + M|(0) where M c (0) is the space-independent 
solution of the gap equation (|J) for a fermi distribution, T represents the temperature, n the chemical 
potential and A the cutoff. The minus sign in the exponentials should be replaced by a plus sign 
in the case of a test anti-particle. In this work we have taken initial quark and energy densities of 
po = B/^tttq = 1.87 fm -3 (which leads to B = 9 for vq = 1.05 fm ), and eo = 1.57 Gev fm~ 3 
respectively. These values correspond, via the definition of quark and energy densities in the presence 
of a cutoff 

2N c N f f A 3 / 1 1 \ 

P °~ (2vr3) J dP \l +e (E P -»)/T l + e {E p+ »)/T) 

2N c N f f A o „ / 1 1 ^ " ' 

eo = "77; t~T" J d 6 p E p I - ; IJ? _, A/rr + 



(2tt 3 ) J 1 p VI + e(^"^/ T 1 + e ^+^l T , 
together with the gap equation (|2|), to the following values for the temperature, the chemical potential 

and the constituent quark mass respectively: T = 240 MeV, fi = 200 MeV and M c (0) = 33.5 MeV. 

The numerical values of A and A for the test-quarks and test anti-quarks respectively have been fixed 

by fitting, at t = 0, to po and eo given above with the result: A = 2630 and A = 630. According to 

Eq. (HJ) the difference A — A determines how many test (anti) quarks one employs for a physical (anti) 



quark [10]. At this stage we want to stress that the well known expressions of the quark and energy 
densities in the case of a fermi distribution with a vanishingly small quark mass: 

P0 = N c N f p(7T 2 + (f3p) 2 )/37T 2 p 2 

e = N c N f (21C(4) + 6(/3/x) 2 C(2) + i(/3/i) 4 )/2vr 2 /3 4 , 
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lead to T = 1//3 = 185 MeV and fi = 190 MeV. The difference between these values and those given 
above is due to the cutoff. To be consistent we have considered in this work the quantities calculated 
with the cutoff A (c.f., Eqs. (||) ). 

For the numerical implementation of our equations, the evolution of the i-th test particle or test 
anti-particle (c.f., Eqs(|5])) is represented by the following algorithm 

Pi (t + St) = Pi (t -St) -2 5t( VM c {r h t)) M c {n,t)/Ei{t) 
r i (t + 5t) = r i (t-5t) + 2 5t Pi (t)/E i (t) 

where 5t = 0.05 fm/c is the mesh in time. At the time t, we determine the constituent mass M c (cell, t) 
(0) at each cell of size a of the space lattice (a being the space mesh-size) according to 

Q ^cell 1 A cell -i 

M c (cell, t) = m + g*M c (cell, t) (^I A (t) ("■) 

where A ce u (A ce u) are all the test particles (anti-particles) which are inside the considered cell, and 
/a = 7rM^(sinh(2a) — 2a) with a = log[ -p- + yj -j^? +T ]. In order to test the accuracy of our numerical 
method, we have checked the energy conservation. We have found that the total energy of the system 
changes by less than 5 percent for 40 iterations in time with a 0.05 fm/c mesh-time. 

Before we come to the numerical results it is instructive to discuss our expectations using a sim- 
plified model. Let us assume for the moment that we have a freely expanding nonrelativistic plasma 
at a temperature T which is initially confined in a sphere in coordinate space. The (Maxwell) velocity 



distribution dN/dv is zero for v = and has a pronounced peak around v = y/2T/m. Without any 
interaction the many particles with about that velocity stay together, i.e. the maximum of the density 
move with v outwards. For radii considerably larger or smaller than r(t) = vt, (t ^> 0) the density is 
lower. Thus the particles are concentrated in a shell. The density in that shell decreases proportional 
to r~ 2 (t). 

Now, what happens if we switch on the interaction. The interaction acts like a repulsive momentum 
dependent potential V = a{p) ■ g(p) where a increases with decreasing density. Since the potential 
is more repulsive in the low density region as compared to the high density region there is a force 
which accelerates the particles towards high density. Thus the high density zone can maintain a 
higher density as compared to the free streaming on expense of the surrounding lower density regions. 
Finally it is energetically favourable to give up the isotropy in <p and 9 and to build local fluctuations 
(droplets) of high density which are separated by regions of lower density. However, this is beyond 
the limits of the predictive power of our one body theory even if we observe these fluctuations due to 
the finite number of test particles. 
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We come now to the numerical results. Fig.l shows the behaviour of the constituent quark mass 
M c as a function of the radius for different times. For high density, i.e. in the plasma phase, M c 
is equal to the bare mass of the quark. At zero density, after the phase transition has taken place, 
M c is equal to the constituent quark mass of ~ 300 MeV. Initially we see that all particles have a 
mass of 33.5 MeV which can be considered as the bare mass. The fact that it is not exactly equal 
to 4 MeV is due to the finite density. Hence the system is in the quark phase. Then the system 
expands and the mass of the particles close to the surface comes close to the constituent quark mass. 
At t = lfm/c we see that the particles have left the center, and therefore the mass has increased also 
there to the constituent quark level. We have now an expanding shell. Before and after the shell 
we have constituent quarks while in the interior we have the plasma. During the expansion the peak 
density decreases and therefore M c increases. The plasma in the interior of the shell starts to make 
the phase transition. This phase transition is almost accomplished after 1.9 fm/c. 

Hence the result of the numerical calculations follows qualitatively the consideration explained 
above. The following scenario of the phase transition emerges. If initially a plasma is created in a 
heavy ion collisions, it is, due to geometry, confined in a limited space region. The high temperature 
of the plasma leads to a fast expansion of the system. The phase transition takes place while the 
system expands. There is a density front traveling outwards. The density is fed by the interaction of 
the quarks which pulls the quarks into the high density region but it also decreases with about r~ 2 (t) 
due to the expansion. Most probably then the system gives up the spherical isotropy and forms quark 
droplets. However, they will not be stable but will hadronize. Obviously this scenario is quite different 
from that expected for a system which keeps global thermal equilibration during the expansion. 

Is this phase transition specific for the NJL Lagrangian or genuine for a phase transition in an 
expanding system? A final answer is not at hand yet but there exists another phase transition in nuclear 



physics which has been investigated in detail recently [11] : The liquid (fragments)-gas (nucleons) 
transition during the expansion of a fireball created at much lower energy ( E^ eam ~ 100 MeV/N) in 
heavy ion reactions. There, for a density smaller than the normal nuclear matter density, we have as 
well the situation that the potential is more attractive at higher density as compared to lower density. 
In the expansion of this system a similar droplet formation has been found in Ref. |ll| if one describes 
the expansion by a 1 - body transport theory as done here. Thus, there is evidence that our results 
are rather general. More realistic iV-body calculations confirm this general structure and show that 
indeed many droplets are formed. 

Calculations including collision terms are under way. Recent calculations of the hadronization 
cross section of a quark plasma fl^] offers as well the possibility to model the chiral phase transition 
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itself. The results will be reported in a forthcoming publication ||. 
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Figure captions 

FIG. 1. The constituent quark mass M c (c.f., Eqs. JTT| ) ) as a function of the radius for 
different times. The initial temperature and density of the system are T = 240 MeV and po = 1-87 
£ m -3 respectively. 
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